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Abstract 
We investigate the connection between energy level crossings in integrable systems and their integrability, 
i.e. the existence of a set of non-trivial integrals of motion. In particular, we consider a general quantum 
Hamiltonian linear in the coupling u, H{u) = T+uV, and require that it have the maximum possible number 
of nontrivial commuting partners also linear in u. We demonstrate how this commutation requirement alone 
leads to: 1) an exact solution for the energy spectrum and 2) level crossings, which are always present 
in these Hamiltonians in violation of the Wigner-von Neumann non-crossing rule. Moreover, we construct 
these Hamiltonians explicitly by resolving the above commutation requirement and show their equivalence 
to a sector of Gaudin magnets (central spin Hamiltonians). By contrast, fewer than the maximum number 
of conservation laws does not guarantee level crossings. 
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INTRODUCTION 



Level crossings - the emergence of degeneracies in a physical system at a certain value of 



some tuned system coupling - underly a myriac 



of compelling phenomena, including anomalies in 



relaxation rates[l|], the onset of quantum chaosj^j, quantum phase transitions [3], Berry's phasey,|5|] 
etc. It is widely believed that these degeneracies can often be understood in terms of a certain 
underlying symmetry. However, in many cases this connection between symmetry and degeneracy 
remains mysterious. This is especially true for quantum integrable systems, e.g. the Id Hubbard, 
anisotropic Heisenberg, reduced BCS models etc. These systems are long known to display an 
abundance of level crossings [6|-l9|, see Fig. [H in violation of the famous Wigner-von Neumann 
non-crossing rulejlCll4l6l| and with no convincing symmetry explanation. 

In this paper we derive the existence of level crossings and an exact solution for a general 
parameter-dependent quantum Hamiltonian from its integrability. Our work has been inspired 
in part by Refs. llT] and |l8| and especially Shastry's paper [19(], which opened up a new, purely 



algebraic perspective on quantum integrable models independent of Bethe's Ansatz. In Hamiltonian 
mechanics the integrability of a system with n degrees of freedom is usually understood as the 
existence of a maximum number (n) of Poisson commuting independent invariants. Then, a well- 
known theorem due to Liouville and Arnold guarantees that the equations of motion can be solved 
by quadratures 



20( 1 . There is no similarly accepted notion of quantum integrability, especially in 
finite dimensional systems, e.g. discrete lattice models in condensed matter physics where the state 
space is generally finite. In particular, it is often unclear what constitutes an independent integral 
and what is the natural notion of the number of degrees of freedom. Nevertheless, it turns out that 
these difficulties can be circumvented if one restricts the manner in which the integrals of motion 
depend on the coupling. 

For concreteness, let us consider Hamiltonians linear in the coupling u. As we are interested 
in discrete energy spectra, we assume that the Hamiltonian can be represented by an A^ x A'^ 
matrix. Following the classical notion of integrability, we require the existence of the maximum 
possible number of independent (see below) mutually commuting integrals, [H^{u),H^ (u)] = 0, 
where H^{u) = T^ + uV^ are Hermitian operators. One of them is the Hamiltonian itself, e.g. 
H^{u) = H{u). Using this commutation requirement alone, we derive an exact solution for the 
spectrum of each H^{u), which can be viewed as an extension of the Liouville- Arnold theorem to 
quantum Hamiltonians. Moreover, we are able to demonstrate that the eigenvalues of W{u) are 
necessarily degenerate at a discrete set of values of u. 



First, we solve the nonlinear commutation relations [H^{u),H^{u)] = to obtain each H^{u) 
explicitly, see below. Interestingly, it tm'ns out that these maximally commuting (or simply max- 
imal) operators H^(u) can be mapped to exactly solvable Gaudin magnets 17|, |2l| (central spin 
Hamiltonians) . The latter describe a localized spin in a magnetic field B = _u_interacting with 



A^ — 1 "environmental" spins and have a variety of physical applications [22l427l|. The mapping 
to Gaudin magnets allows us to obtain the exact solution for the eigenvalues and eigenfunc- 
tions of W{u) and analyze its level crossings. The total number of pairwise crossings varies as 
Mc = {N — 1){N — 2)/2 — 2K > 0, where i^ is a positive integer such that Mc > 1. For instance, 
A^ = 5 maximal operators display 2, 4, and 6 instances of level crossings. On the other hand, we 
demonstrate that Hamiltonians having less than the maximum number of n-dependent integrals 



can have no level crossings. Further, we show that the coup 



ing-dependent commuting matrices 



obtained within the approach developed by Shastry in Ref. [191 (where it was also conjectured that 
these matrices always have crossings, see also Ref. .18*) are maximal, even though our constructions 
are quite different. 

Pairwise crossings of energy levels are usually understood in the context of the Wigner-von 
Neumann non-crossing rule. This rule initially suggested by Hund[lO| and justified by Wigner and 



von Neumann [ill] has thereafter seen restatements and refinements by a number of authors|l2l4l6|. 
It states that eigenstates of the same symmetry do not cross as a function of a single coupling 
parameter. This can be seen, for example, from the following argument. Suppose two energy 
levels El (n) and £'2 (u) of H (u) come close at a certain u = uq. Expanding in a vicinity of uq: 



H (u) ^ H (uq) + (n — Uq) V (uq) and using ordinary perturbation theory, we obtain [28|] 

where A (uq) = Ei (uq) — E2 (no), V12 (uq) is the matrix element of the perturbation V (uq) between 
states |1) and |2), and F (uq) represents the contribution of the remaining states. We see from 
Eq. ([1]) that as the two levels approach, A (uq) -^ 0, infinite repulsion sets in, preventing them 
from crossing. This is indeed what takes place in the absence of symmetry - energy levels repel 
see Fig. [2j The situation changes if the Hamiltonian H (n) possesses a u-independent symmetry 
S, i.e. [H {u) , S] =0. This can be a spatial rotation, translational invariance, internal space 
reconfiguration, etc. Because S does not depend on the coupling u, it commutes with H (uq) and 
V {uq) individually. Evaluating the matrix element of [V{uq),S] between states |1) and |2), we 
obtain ^12(^0) = for any uq as long as |1) and |2) have different symmetry, si 7^ S2, where 
5|1) = si |1) and S\2) = §2 1 2). Thus, while levels of different symmetry can cross, crossings of 
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FIG. 1: Energies of Id Hubbard model on six cites characterized by a complete set of quantum numbers, i.e. 
all levels have the same u-independent symmetry, see Refs. |6| and [iSl The energies are in units oi U — 4r 
plotted as functions of m = U/{U — 4r), where U is the strength of the Coulomb repulsion and T < is the 
hopping matrix element. The parameter u varies from to 1 as [/ goes from to oo. Note that, in violation 
of the Wigner-von Neumann non-crossing rule, we see a profusion of level crossings for states of the same 
symmetry. 



levels of the same symmetry are prohibited. 

Unfortunately, this basic argument does not extend to quantum integrable Hamiltonians H (u), 
which typically violate the non-crossing rule. Indeed, these systems show crossings of energy lev- 
els that have the same quantum numbers for all u-independent symmetry [a-lSl 5 see e.g. Fig. [TJ 
Integrable Hamiltonians are known to have special coupling dependent conserved currents, "dy- 
namical symmetries", in addition to u-independent symmetries. It is tempting to attribute these 
crossings to such symmetries. On the other hand, it is crucial for the validity of the non-crossing 
rule that the symmetry S be u- independent. Indeed, consider an integrable Hamiltonian acting 
on a finite-dimensional space, e.g. a lattice model with a finite number of sites. Let H (u) be one 
of its blocks characterized by the same quantum numbers for a complete set of mutually commut- 
ing M-independent symmetries and let H (u) be the corresponding block of one of the conserved 
currents 



[H{u),H{u)]=0 for all u. 



(2) 



Due to the u-dependence, H (uq) no longer commutes with V (uq) separately and, therefore, the 
above argument lifting the level repulsion does not hold. At the same time, given a crossing 
one can always artificially engineer a "conserved current" that commutes with H{u). Therefore, 




1.5L 



FIG. 2: Energy levels of a 14 x 14 Hamiltonian H (u) = A + uB, where independent matrix elements of A 
and B are uniformly distributed random numbers. Note that though levels do approach one another closely, 
they never cross. A typical Hamiltonian with no u-independent symmetry exhibits such level repulsion, see 
Eq. (IT|), in contradistinction to what takes place in integrable systems, see Fig. [TJ Here and throughout this 
paper we use Mathematica program to perform numerical tests and plot the results. 



restrictions on the form of H (u) are necessary to make meaningful contact between symmetries 
and degeneracies. 

To be specific, let H{u) = T + uV and H{u) = T + uV be Hermitian operators acting on an 
A^-dimensional space, i.e. they can be represented by A^ x A^ matrices. Eq. ([2|) implies 



[T,T] = [V,V] = 0, [T,V] = [T,V]. 



(3) 



For any linear H (u) there are always trivially related commuting partners H (u) = aH (n) + 
(6 + cu)/, where / is an identity matrix. However, the requirement that Eq. ([3]) have nontrivial 
solutions leads to a set of nonlinear constraints that severely restrict the matrix elements of both 
H{u) and H{u). For example, for A^ = 3 eliminating T and V from Eq. ([3]), one obtains a single 
nonlinear constraint on the matrix elements of H{u) [l8|]. In view of the preceding discussion 
regarding the prevalence of level crossings in integrable models, a natural question is whether these 
constraints, i.e. the existence of a nontrivial H{u), imply crossings in the spectrum of H{u) and 
vice versa. This is indeed the case for A^ = 3. Specifically, one can show that 3x3 matrices 
H{u) = T + uV that have nontrivial commuting partners also have a level crossing and vice 
versa[18:]. However, this is no longer true for A^ > 4 - Eq. ([3]) does not necessarily lead to level 
crossings. Moreover, crossings occur even in the absence of nontrivial partners and u-independent 
symmetries, see below. We see that a single dynamical symmetry is insufficient to explain level 



crossings. On the other hand, quantum integrable Haniiltonians typically have more than one 
coupling dependent commuting operator. In fact, as we show below, the maximum possible number 
- which turns out to be A^ - of integrals are necessary to ensure level crossings. 

We define the set of maximally commuting Hamiltonians as a vector space, Ai, formed by 
A^ > 3 Hermitian, mutually commuting N x N matrices W (u) = T^ + uV^ together with the 
N X N identity matrix /, 

[W (n) , H^ {u)] = 0, for all u and ij = I, . . . , N, (4) 

where u is a real parameter. Operators W{u) are assumed to be independent in that matrices V^ 
are linearly independent, i.e. X]j=i ^i^* = iff Ci = for all i (equivalently one can require that 
T* be linearly independent). In addition, W (n) are taken to have no li-independent symmetry 



common to all H^ (n)J3[l|. 



t^^al such that [17, H\u)\ = for all -u and i. (5) 

Therefore, an arbitrary element H[u) = T + uV of the vector space M has the form 

N 

H{u) = Y^ diW {u) + al, (6) 

i=l 

where di and a are real numbers. The addition of multiples of the identity affects neither commu- 
tation relations nor level crossings of H{u) and we will often omit the term al in Eq. ([6]). Note 
also that Eq. ([6]) implies that operators H^{u) together with / provide a basis in the vector space 
^A of maximal Hamiltonians. 

The set M is maximal in the sense that any Hermitian H (u) = T + uV that commutes with 
all H^{u) can be written in the form ([6]). Indeed, since V and all V^ mutually commute, see 
Eq. dS]), we can go to their common eigenbasis. In this basis, the N diagonal matrices V^ are N 
linearly independent A^-dimensional vectors and, therefore, there exist real numbers dj such that 
V = J2i=i diV^- The matrix H{u)—J2i=i ^iH^i^) is ^-independent and, since it also commutes with 
all H^{u), it must be of the form al according to Eq. ([5]). Thus, H{u) is of the form ([6]). By a similar 
argument one can show that one of the basic matrices H^{u) can be chosen as H^{u) = (a + ub)I 
with real coefficients a and b. We see that there are A^ — 1 nontrivial independent commuting 
operators. Therefore, the first nontrivial dimensionality is A^ = 3. 

In what follows we begin with the explicit construction of a general, maximally commuting 
Hamiltonian H{u). This is done in Sec. lUby choosing a convenient basis in the vector space M 
and solving Eq. (JH). In Sec. lIIII we establish some useful algebraic properties of H{u). Interestingly, 



it turns out that the product of any two maximally commuting Hamiltonians can be written as a 
linear superposition of such Hamiltonians, i.e. the set Ai has a certain quasi-ring structure. 

^y_commuting Hamiltonians makes it transparent that they 



17 



21[, see Sec. IIVI The latter are A^ quantum spin Hamilto- 



Our parametrization of the maximal 
are related to the Gaudin magnets 
nians 

N 2, 2. 

W{B) = -Bst + Y,'^^^^, z = l,...,iV, (7) 

where the prime indicates that the summation is over k ^ i, B is the z-aligned magnetic field, Sj is 
an operator of spin of magnitude Si, and l/(ei — £k) is the coupling between spins Sj and s^. The 
Hamiltonians /i* form a mutually commuting family 

h\B),h^{B)\=0 for any 5. (8) 

Note also that each h^{B) is invariant under rotations around the z axis, which means the z 
component of the total spin J^ = Yli=i ^l i^ conserved 



h\r 



0, i = l,...,iV. (9) 



As we will see in Sec. IIVI the maximally commuting Hamiltonians ([6]) correspond to the sector of 
Gaudin magnets with J^ = J^ax ~ 1' where J^ax = Z]i=i ^-i i^ ^^^ maximum eigenvalue of J^ . 

In Sec. |V] we employ the mapping to the Gaudin magnets to obtain the exact solution for the 
spectra of maximally commuting Hamiltonians. Using this solution, we analyze the asymptotic 
behavior of the eigenstates in the limits u — > itcxD in Sec. IVI[ Matching the two limits allows us to 
establish the presence of level crossings and to count them. Next, we consider Hamiltonians having 
less than the maximum number of commuting partners. In Sec. IVIIi we construct a family of 4 x 4 
Hermitian operators linear in the coupling u that have a single nontrivial partner as opposed to 
two partners for the A^ = 4 maximal set. It turns out that these submaximal Hamiltonians often 
display no degeneracies at all. Finally, in Appendix we review Shastry's approach to constructing 
commuting real symmetric operators and show that the resulting operators are always maximal. 

II. THE PARAMETRIZATION OF MAXIMALLY COMMUTING HAMILTONIANS 

We begin our analysis by choosing a convenient basis in the vector space of maximally com- 
muting Hamiltonians, which allows us to solve Eq. ([4]) explicitly. The solution yields a convenient 
parametrization for a general maximal Hamiltonian ([6]) and elucidates the algebraic structure of 



these operators. It also makes transparent the relationship between these operators and Gaudin 
magnets ([7]). 

Consider the maximal operators W{u) = T* + uV"^ defined in Eqs. (JH) and ([5]). It follows from 
Eq. ([1]) that all V^ commute with each other, see Eqs. ([3]) and ([2]). By a u-independent unitary 
transformation we go to a basis where all V"^ are diagonal. Since V^ are also linearly independent, 
we can bring them to the following "canonical" form by a linear transformation 

D' = Y,v':v\ (10) 

i 

where v^ are real numbers, D are diagonal with a single nonzero matrix element [D ]jj = D^ = 6jk- 
Next, we introduce a "canonical" basis in the space of maximally commuting operators 

N 

h\u)=E' + uD' = ^v)W{u), i = l,...,N. (11) 

The operators h^{u) have all the properties of maximally commuting Hamiltonians defined in 
Eqs. (jl]) and ([5|) as long as H^{u) do. In particular, 

[h'{u),h^u)]=0, ioralluandi,j = l,...,N. (12) 

It follows from Eq. ^ that a general maximally commuting operator can be written as 

N 

H{u)=T + uV = Y^ dkh!" (n) + al, (13) 



fe=i 



where d^ and a are real numbers. Note that with our choice of D , d^ are the eigenvalues of V . 

To determine H{u) explicitly, we need to solve Eq. ()12p . In terms of D* and E"^ these equations 
read 

[D\D^] =0, [D\E^] = [D\E^ , [E\E^] =0. (14) 

The first equation holds since D^ are diagonal. The second equation in terms of matrix elements is 

(' n* _ n* "i FJ = (' ni _ dj \ p* (-\ k\ 

where Emn is the mn}'^ matrix element of E^ . By construction, the only nonzero matrix element 
of D' is Dl = 1. We see that Emn = as long as m ^ n and m and n do not equal j. Thus, matrix 
E^ only has nonzero elements of the form i?^^ = [E^-) and Emm, where z* denotes the complex 
conjugate of z. Note also by setting m = i and n = j in Eq. (fT5|) that Ej- = —E^- for i / j. 



It remains to solve the last equation in (J14p . Using the above properties of matrix elements of 
E^, we rewrite this equation as 

rpm Tpm 
^11 ^mm jpra Tpm^ i -f- in, n, 

rpm rpm 

jpm Tpm 
jpn _ rpn _ rpm _ rpm , S^ mj jm 



nn -^mm -^mm -^nn ' / _, rpm jpn 

j^m,n "" ~ iJ 



By direct computation, one can show that the following ansatz satisfies Eqs. ()16p : 



ipm _ „t{em-e,) 7m 7i 



7^ 
-im, im 



Ef^ = --^ + ^-, z/m, (17) 



Tpm _ \ ^ H L „ir, 

^mm - - l^ r -c ^ ^ 

., ^m tj 
t^m 



2 



where jj ^ 0, £j, 9m, and ip"^ are real parameters. A nonzero value of ip"^ corresponds to an overall 
shift of the diagonal of £"", which yields a (nonessential) contribution ^"^I to h"^{u). Note that 
7j = also satisfies Eqs. ([16]) but generates matrices with block diagonal structure and, therefore, 
ti-independent symmetry. 

Furthermore, any solution of Eqs. (J16p admits parametrization (J17p . To establish this, it is 
sufficient to show that any choice of the 3A^ — 2 matrix elements E^^ = (-E^) and ii^^ for a 
certain m compatible with Eqs. (|16p corresponds to a set of 3A^ + 1 real parameters, 9j, jj, £j, and 
ip"^. Then, Eqs. (116p ensure that all E^ are of the form (I17p . The extra three parameters are an 
overall scale for £j and 7j, a shift e^ — > e^ + const, and a shift 9j — > 9j + const, which do not affect 
Eqs. (fT7|) . To see the correspondence between the E^^ = (.E^) and £'^ and the 9j,jj,ej, and 
V^™, note that Eqs. (^ yield 



ipm Tpm 

E rnn nm, jpm, /m 
p?ni _ ,/,r?i ~ rnm 'P ' 



Tpm ,/,n 



which, for given E^n = {^rvm) ^^'^ ^mi.^ '^^^ be solved for -0™. In seeking a common denominator, 
we see that it constitutes an N^^ order polynomial 

n {Ef, -i^n-Y. n ^-n^r-(^i; - ^n = 0' (18) 

j n^mj^m,n 

10 



yielding A'^ solutions ip"^. By considering the form of matrix E"^, i.e. that all matrix elements are 
zero save for a row, the corresponding column, and the diagonal, we find that the left hand side of 
Eq. (jlSp is the characteristic polynomial given by det {E^ — ip^I). Consequently, all N solutions 
of Eq. (jlSp are guaranteed to be real as they are the eigenvalues of an explicitly Hermitian matrix. 
Once ip'^ is determined, we can calculate ratios 

7j - e _ J ^ m, 

jj ^ 

where e^*(^^~^™) = E^/E^-, 9m is arbitrary and by a choice of an overall scale we set 7™ = 1- 
Lastly, letting e^ = (by shifting Sj), we have 

1 



e 



1 Ej^.-iprn- 

Each of the A^ solutions ■i/'™' to Eq. (|18p will yield a distinct set {jjjCj}, but by construction (see 
Eq. (fT8|) ) corresponds to the same set {Em^,El^m^Ejj}- 

Now, consider i?™ as defined in Eq. (J17p . The matrix £"" with complex matrix elements 
E^j = (EJ^)*, E^, and E^^ is conjugate to a matrix E"^, i.e. 






where S is a diagonal matrix with entries T^jj = e*^^ = e"^"'.^, {E^^jY = E^jEJ^, E^Pj = E^j'^, 
and 9m is an arbitrary real number. Given a Hermitian E"^, we find that E^™" is necessarily real 
symmetric matrix and S is a conjugating diagonal matrix whose matrix elements are complex 
phases. Thus, all Hermitian maximally commuting operators are matrix conjugate to some real 
symmetric such. Operator spectra are invariant under matrix conjugation and, therefore, it is 
convenient to henceforth limit our discussion to maximally commuting real symmetric matrices, 
and we do so without loss of generality. 

Moreover, as noted belowEq. (J17p . nonzero -0* contributes only a multiple of the identity, ip'^I to 
each /i* (u) , which affects neither the commutation relations nor the level crossings of /i* (n) and their 
linear combinations. Henceforth, we adopt a convenient "gauge choice" V* = for alH = 1, . . . , A'^ 
and 9m = for all m = 1, . . . , A^. With this choice of ip^ and 9m, we derive from Eqs. (|ll|17p and 



11 



the definition of D* the nonzero matrix elements of basic operators h'^{u) 



f 



[/.^(n)],, = -^, ./,, (19) 



[hHu)i,=u-Y: 



7. 



2 



3 



3+1- 

Note that 

AT 

^ U [u) = ul. (20) 

Expressions (J19p for matrix elements constitute a complete, explicit solution of commutation re- 
lations (J12p or equivalently (j3|) for maximally commuting Hamiltonians. Different choices of pa- 
rameters 7j and Ej (factoring out overall scale of jj and Ej , a total shift of all Ej , and the "gauge 
freedom" discussed above) yield distinct families of such Hamiltonians. 



Eq. ()19p also determines matrix elements 



29I ] of a general maximally commuting operator (fT3 



-m c^n , 



(21) 

j^m ^^"^ ^^^ 

Let us also note that a convenient approach to producing nontrivial solutions of Eq. ([2]) was 
developed by Shastry in Ref. [l9|. Interestingly, these solutions turn out to be essentially equivalent 



to the maximally commuting set constructed in this section, see Appendix for details. 

III. ALGEBRAIC PROPERTIES OF MAXIMAL HAMILTONIANS AND AN UPPER 
BOUND ON THE NUMBER OF LEVEL CROSSINGS 

The above parametrization makes transparent a beautiful property of maximal Hamiltonians - 
the product of two maximal operators is itself the u-dependent sum of maximal operators. This 
property, as we demonstrate in this section, allows one to express a general maximal Hamiltonian 
H(u) as a polynomial in another such Hamiltonian H{u). We employ this polynomial expansion 
to determine the maximum number of level crossings in the eigenvalue spectrum of H{u). 



12 



First, we express the product of two basic maximally commuting operators h^(u) and h^ {u) in 
terms of u-dependent linear combinations of h^{u). Using Eq. (|19p . one can show that 



h' (u) ■ y (n) 



72 ^2 

^ -h'{u)+ '' 



En 



N 



h} (u) . h^ (U) = ^ 



fc=l 



E 



En 



Tn 



-h^{u), i^j, 



(22) 



Efc 



h''{u) 



Now consider two general maximally commuting Hamiltonians (jl3|) 

N N 

H{u) = Y,dkh''{u), H{u)=Y,dkh'{u), 
fe=i fc=i 

where without loss of generality we dropped multiples of identity in Eq. p3|) . From Eq. 

derive 



TV 



H{u)-Hiu) = Y, 



k=l 



U 



dkdk — / ^ 

rriy^k 



1m (dk - dm) {dk - dr, 

Ek £'m 



h'' (n) 



(23) 



we 



(24) 



This quasi-ring structure - so called because, while the sum of maximal Hamiltonians is maximal, 
the product is a u-dependent sum of such and, therefore, not generally linear in u and not strictly 
a maximal operator ~ suggests a means of representing an element of a commuting maximal family 
by any other, see Eq. (j25p below. 

A typical maximal Hamiltonian H{u) can be degenerate only at discrete values of u. Note 
that the only alternative to the discrete (possibly empty) set is a permanent degeneracy - when 
two eigenvalues of H(u) coincide at all u [31,]. Permanent degeneracies do not occur for a generic 
choice of df^ in Eq. (j23|) . Indeed, recall that d^ are the eigenvalues of V (see below Eq. (|13p ). Since 
the eigenvalues of H(u) = T + uV tend to those of uV for large u, the spectrum of H{u) is not 
degenerate as long as the dk are distinct and V is itself nondegenerate. 

Consider H{u) at any u where it is nondegenerate. Any element of its commutant - the set of 
all real symmetric operators that commute with H{u) - can be expressed as a polynomial in H{u) 
of the order A^ — 1, i.e. 

N-l 

H{u) = J2Pa{u)H"{u), (25) 

where H" (u) = [H (u)]" and, as we will see shortly. Pa {u) are rational functions of u. To see 
that H{u) can be indeed written in terms of powers of H{u), consider Eq. (j25p in the common 
eigenbasis of commuting operators H{u) and H{u) at a given u. Since eigenvalues ujm of H{u) 



13 



are A^ distinct real numbers, one can always find a polynomial Ri\f-i{uj) = X^^Jg ^a'^" of order 
A^ — 1 with A^ real coefficients Pa so that Ris[-i{uJm) = ^mi where ujm are the eigenvalues of H[u). 
Indeed, the equations RN-i{oJm) = ^m are linear in Pa with a nonzero determinant. 
Next, we observe from Eqs. ([23|) and (p^ that 

N 

H-{u) = Y,Qtiu)h'{u), (26) 

fc=i 

where for a > 1 Q" (u) is an a — 1 order polynomial in u determined by recursively applying 
Eq. (f24l) and Ql = 1/u as follows from Eq. ([20]) . Plugging Eq. ([26]) into Eq. (|25]) and using the 
second equation in (p3|) . we obtain 

N N~l N 

^ Y, Pa{u) Ql{u) h\u) =Y,dk h\u). (27) 

k=l 0=0 fc=l 

Since h'^ {u) are linearly independent at any it / 0, i.e. J2k=i fki"^) ^^i''^) = if and only if 



fk{u) = for all -u 7^ 



32], Eq. (HZD becomes 

Af-l _ 

J^ P„ (n) QH^^) = 4, k = l,...,N. (28) 

0=0 

Note that because Q'^ (u) are rational functions in n, PQ,(n) are also rational functions. 

Because H{u) is arbitrary, Eqs. ()28p should have solutions for P^ (n) for any dk as long as H{u) 
is nondegenerate. On the other hand, solutions cease to exist if and only if det[(5^(u)] = 0, where 
Q^{u) is regarded as the ak matrix element of an A^ x A^ matrix. Using the fact that Q^ (n) is a 
polynomial in u of degree a — 1 for a > 1 and Q\ = 1/u, one can show that det[(5^(n)] = V{u)/u, 
where P(n) is a polynomial in u of order "^^Z^m. The real roots of the equation V{u^) = are 
the values of n = {u^} where H{u) is degenerate. Thus, the maximum possible number of level 
crossings in the eigenvalue spectrum of an A^ x A^ maximally commuting Hamiltonian is 

M»" = '^-'>f-^'. (29) 

The polynomial V{u) is of real coefficients and, therefore, its complex roots come in conjugate pairs. 
Consequently, the number of real roots of V{u) falls from the maximum M^^^ in decrements of 
two. This enforces a parity such that the number of real roots is odd for integers of the form 4?7i, 
4m + 1 and even for integers 4?ti- + 2, 4?7i + 3, m G N. Ostensibly, when real roots of V{u) are 
degenerate their number need not correspond to the number of distinct crossings. In principle, a 
multiple real root oiV{u) could correspond to a single pairwise crossing. Numerically, however, we 
have observed that such multiplicities occur only when more than two levels cross simultaneously, 
i.e. at the same value of u. 

14 



IV. MAPPING TO THE GAUDIN MAGNETS 

In this section, we show that maximally commuting Hamiltonians h^{u) are equivalent to the 
Gaudin magnets, 



h'{B) 



-Bs 



N 2. 2. 

z , sr^' Si- Sk 



+ E 



,N, 



k=l 



Sk 



(30) 



in the next to highest weight sector, J^ = J„ 



1, where J^ is the z projection of the total spin, 



Si = {sf,'sf,'S|} is the i^^ spin 3-vector of magnitude Si and [sf^sA = eajs-ysjdij. This ma pping is 



very usefu. 



M, 



35M 



as Gaudin magnets (central spin Hamiltonians) have been extensively studied 



iU- 



17| 



For example, an exact solution for the eigenstates and eigenvalues is available 



21 



13, 






We employ it in subsequent sections to obtain the spectra of maximally commuting Hamiltonians 
and to analyze their level crossings. This mapping also implies that all our conclusions regarding 
maximal Hamiltonians, e.g. the presence and the number of level crossings, quasi-ring structure 
(j22p etc. can be immediately transferred to the cor resp on ding sector of Gaudin magnets and their 



3S 



40l |. At the same time, other sectors of 



derivative models, such as the reduced BCS mode! 

the Gaudin model as well as more general models [36| of which it is a particular case can provide 

examples of Hamiltonians with less then the maximum number of commuting partners. 

Since Gaudin magnets ([30]) commute with the z projection of the total spin J^, see Eq. ([9]), they 
are block-diagonal in any basis where j^ has a definite value. Different blocks can be labeled by 
the eigenvalues of J^. Consider the sector J^ = J^ax ~ 1, where J^ax = Yli=i Si is the maximum 
eigenvalue of J^. It is populated by A^ basic states 



/2il 



k 



iV, 



(31) 



where |0) is the highest weight state J^ = J^i^^, i.e. s^|0) = for all k, and the highest weight Sk for 
each spin s^ is given by s||0) = s^jO). Therefore, Gaudin Hamiltonians (j30|) are N commuting real 
symmetric N x N matrices in this sector. Since there is also no obvious i?-independent symmetry 
(J^ oc / within a given sector), the h^{B) appear to be good candidates for a maximally commuting 
set. 

To check this, let us evaluate the nonvanishing matrix elements of h}[B) given by Eq. (j30p in 
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the normalized basis (1311). We obtain 



{i\h\B)\j) = ^^^, j^i, 



£i Ej 



m {B) ij> 



Ei Sj 



+ 



-Bsi + Y, 



k^i 



£i — £k 



J /«, 



(32) 



{i\h\B)\^ = B-Y,-^ + 



k^i 



. ^i — £k 



-Bsi + Y, 



k^i 



£i — £k 



Comparing these expressions to matrix elements of h^{u) in Eq. ()19|) . we observe that with the 
identifications B = u and Sk = 7I the two matrices differ only by a multiple of an identity matrix 
ip^I, where 



^' 



-Bsi + Y, 



k^i 



^i ^k 



(33) 



Recall that we arbitrarily selected a "gauge" -0* = for maximally commuting Hamiltonians h^{u), 
see Eq. (J17p and the text above Eq. (J19p . This constant overall shift of all eigenvalues of /i*(u) 
affects neither its eigenstates nor the degeneracies. 

Thus, we see that Gaudin Hamiltonians (j30p in the next to highest weight sector J^ = J^a,x ~ ^ 
are equivalent to basic maximal Hamiltonians h^{u) with 



u = B, 7I = Sk, 



(34) 



and vice versa. Note that the magnitudes of quantum spins, s^, take half-integer values for finite 
dimensional representations of the spin su (2) algebras, while 7^ are arbitrary real numbers. We 
believe that this restriction can be lifted by moving to an appropriate infinite dimensional repre- 
sentations of the su (2)s, where the highest weight states are still well defined but Sk take arbitrary 
real values 4]|. Indeed, we have verified that, at least in our sector J^ = J^^^^ — 1, in all expressions 
for the eigenvalues and eigenstates of h'^{B) (see below) the replacements B ^ u and s^ -^ 7^ with 
arbitrary real 7^ produce the correct corresponding eigenvalues and eigenstates of h^{u). 



V. EXACT SOLUTION FOR THE SPECTRA OF MAXIMAL HAMILTONIANS 

A particularly useful consequence of the mapping ()34p between Gaudin magnets h^(B) and 
maximally commuting Hamiltonians /i*(ti) is that one can obtain the exact solution for h'^{u) by 
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importing the known exact solution for the spectra of h'^{B) jlTI . |2l|. The latter has been derived 



both from the properties of the Gaudin algebra|l7l| and by Bethe's Ansatz 36l |. 

The exact eigenvalues of the Gaudin Hamiltonian (130p . h^{B), in the next to highest weight 



sector J^ 



•^max ~ 1 ^^^ 



[Kn) 



G 



x: 



G 



Ei 



+ V' 



(35) 



where V* is the overall shift of all eigenvalues given by Eq. (j33p and x^ are the solutions of the 
following equation: 

TV 



B 



E 

fc=i 



„G 



£k 



(36) 



Note that if this equation is brought to the common denominator, the numerator becomes a 
polynomial of order A^ in x^. Therefore, there are A^ solutions for x^ and A^ eigenvalues (I35p as 
it should be since there are A^ states in this sector, see Eq. (I3ip . The unnormalized eigenstates 
(common to all h^{B)) corresponding to eigenvalues (j35p are 

f:#^. (37) 



A^ 



'G 



x: 



£k 



k=i "* 
where the basic states \k) have been introduced in Eq. (j3ip . A concise derivation of Eqs. (J35l40p 



and (j37p can be found in Refs. 
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21 



, and 



37 



Using the mapping (j34p between basic maximal operators h'^{u) and Gaudin Hamiltonians, we 
obtain from Eq. (I35p the energies of /i*(n) 

v2 



A^ 



X 



(38) 



Note that we set the overall shift -i/;* = in accordance to the discussion surrounding Eq. (133p . The 
corresponding common eigenstates of all h^{u) are 

N 



|A^) = Y. 



fc=l 






where |A;) now stands for a basic vector for matrices K'{u), i.e. its j component is \k)j 
Eqs. ([38]) and (f371) Xm are solutions of the following equation: 



N 

E 

fc=i 



ll 



£k 



= fiXr. 



m 



l,...,N, 



(39) 

Sjk- In 

(40) 



which follows from Eq. (j36l) . That Eqs. (I38|39p and (I40p yield the correct spectrum of h'^{u) can be 
verified directly using the matrix form ()19p of /i* (u) . Finally, using Eq. (|23p , we derive the energies 
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Il 



X2 



X3 




X4^ 




FIG. 3: A plot of f{x) — J2k=i Tfel-^ ~ £fc)~^ for A^ = 6. Solutions of the equation f{xm) = u determine the 
energies (j4l]) and eigenstates (|39|) of a general maximal Hamiltonian H{u). There are N points, a;,„ with 
m = I, . . . ,N, where y = f{x) intersects y = u (dashed horizontal line) yielding A^ = 6 eigenstates. Note 
that Sm < a^m < £m+i except in the case of xn for which we have Sn < xn for u > and xn < £n+i = £i 
for u < 0. Furthermore, we see that Xm — >■ £m as w — > oo and Xm —^ £m+i as m — > — oo. This observation 
allows us to determine the behavior of the energies in the u — > ±cx) limits, see Eq. (|44p . 



of a general maximally commuting Hamiltonian, H{u) = T + uV, 



N 



Wn 






k=l 



£k 



(41) 



The corresponding eigenstates are still given by Eq. (|39p . 

Let us analyze the flow of eigenvalues tOm of H{u) with u and determine their behavior in the 
u -^ iboo limits. Consider Eq. (00]). Let e^ be ordered as ei < £2 < • • • < ^N- The right hand side 

k+i- 



of Eq. (jlO]) is plotted in Fig. [3l Note that f{x) — > +00 as x ^ e^ and f{x) — > —00 as x — > e 



It follows that the equation u = f{xm) has a real root between e^ and £^+1 for any k. Let us 
number the roots Xm so that Em < Xm < ^m+i- Note from Fig. [3] that for the last root X]y we have 
£m < xn for u > and xn < £n+i = £1 for u < 0, where from now on we identify indices m and 
m + N that differ by a multiple of N . Further, observe from Fig. [3] that Xm — > em as u — > +00. 
In this limit k = m term dominates Eqs. (I40p and (I4ip and we obtain Jm/{xm — ^m) ~ u and 
Wm ~ dm7m/(a;m " ^m) ~ "U-dm- Similarly, for u -^ -cx) we have Xm -^ £m+i and uJm ~ "Ufim+i- 
Therefore, 



w^ 



|u|dm+i as tt — > —00, 



Wr, 



|li|dm as li ^ +CXD. 



(42) 
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At this point it is convenient to rescale the Haniiltonian 

H',, .^ L^ ,43, 

VM + 1 V u + 1 

Note that this does not affect the level crossings, i.e. H'{u) and H{u) have crossings (if any) at 
the same values of u. Eq. (j42l) implies 



u^ —>■ —dm+i as u — > — oo, o;^ -^ dm as u ^ +oo, (44) 

where u;[„ is the eigenvalue of H'{u) corresponding to the eigenstate |Am)- Recall that d^ are 
the eigenvalues of F, see the text below Eq. (fT3]) . We see from Eq. (|13|) that the eigenvalues of 
H'{u) indeed should tend to dfc in u ^ ±00 limits consistent with Eq. (j44p . The latter equation 
however provides much more detailed information - it shows to which particular dk the eigenvalue 
corresponding to a given eigenvector tends in each limit. We will use Eq. (|44p in the next section 
to study the crossings of energy levels of a general maximally commuting Haniiltonian H{u). 

VI. LEVEL CROSSINGS 

In this section, we establish the presence of energy level crossings in the spectrum of an arbitrary 
maximally commuting Hamiltonian H{u) = T + uV ([6|). This provides an explanation of the level 
crossing phenomenon in the absence of any n-independent symmetry based solely on the fact that 
H{u) has the maximum possible number of independent commuting partners, see the text above 
Eq. ^. Further, we determine the number of level crossings as it depends on the ordering of the 
eigenvalues dk of the perturbation operator V and argue that this number takes values 

Me = ^^ ~ ^^^^ ~ ^^ - 2K, K = 0,1,...,K^^, (45) 

where A^ is the dimensionality of the state space of H{u) and ATmax is the integer part (floor) of 
{N — 1){N — 2)/4 — 1/2. For example, A^ = 3 maximally commuting operators have a single level 
crossing, while for A^ = 6 we have Mc = 2, 4, 6, 8, and 10. We also develop an approach that allows 
us to readily predict the minimum allowed number of crossings for a given H(u) from the ordering 
of dk- 

Consider a Hamiltonian (not necessarily belonging to any commuting family) that depends on 
a real parameter u. Suppose |nj) are its eigenstates and E^_^ < E^,^ < . . . are the corresponding 
energies at large negative u. There is only one way to avoid crossings - the order of eigenvalues 
E^_ at u — > 00 must be exactly the same as that at li — > —00, i.e. E^_^ < E^^ < . . . This is what 
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happens with a typical Hamiltonian in agreement with the Wigner-von Neumann non-crossing 
rule, Fig. [5J If, on the other hand, the relative order of any two energies changes, at least one 
level crossing must occur. For example, E~_^ < E~^ and E^^ > E^^ means that the difference 
Eniiu) — En4^{u) changes sign as u evolves from — oo to co. By continuity this implies a crossing of 
levels corresponding to eigenstates |ni) and |n4) at a certain value of u. This is observed in blocks 
of quantum integrable Hamiltonians characterized by the same u-independent symmetry, see e.g. 
Fig. [TJ Numerical spectra of maximal Hamiltonians display the same behavior. Figs. [5] and [8^. 

Now let us turn our attention to an arbitrary maximally commuting Hamiltonian H{u). In the 
previous section we have established the behavior of its energies in u ^ ±00 limits. It follows from 
Eq. (j44p that the energy level of H{u) = T + uV (with appropriate rescaling (j43p ) that starts from 
—dk at It — > —00 ends at dk-i at ti — > 00. Symbolically, this can be represented by 



k 



k-1, (modiV). 



(46) 



Note that we cannot fix an ordering of dk without loss of generality, as the d^ correspond to Sk, 
see e.g. Eq. (jlT]) . and we have already fixed the order of Ek so that ei < 82 < ■ ■ ■ < £n- First, we 
assume that all dk are distinct as is generally the case. Eq. (j46p implies that the flow of energy 
levels from u = — cxd to u = cxd can be schematically depicted using the following rules: 



1. Create two columns in which {—dk} and {dk} are both in descending order and replace each 
dk with its lower index k, i.e. 



-di 
-dj 

-di 



di 



I 
3 

I 
m 



m 
I 

3 
i 



(47) 



2. Draw a line connecting j in the left column to j — 1 in the right. These lines represent energy 
levels of H[u). Consequently, their crossings imply crossings of the corresponding energy 
levels of H{u). 

An example of an energy level diagram generated using the above prescription for A^ = 5 is shown 
in Fig. m It corresponds to the ordering ^4 > ds > ^3 > di > ^2 and predicts six level crossings. It 
also specifies which levels cross, e.g. the top level connecting 2 and 1 crosses with the next in energy 
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FIG. 4: Schematic energy level diagram for general 5x5 maximal Hamiltonians H{u) = T + uV drawn 
using the rules (|47|) . The eigenvalues dk of V are ordered such that d^ > d^ > ds > di > d2. The diagram 
shows six level crossings for this ordering and specifies which levels cross, e.g. the level 4^3 first crosses 
3 — > 2 and then 2^1. Compare to Fig. [5] and note that the crossings predicted by the above diagram are 
exactly the same as actual numerical crossings for this ordering. Note also that according to Eq. (|15|) six 
is the maximum allowed number of crossings for this ordering of d^. and multiple crossings of the same two 
levels are therefore forbidden for this ordering. 



level connecting 1 and 5. We see that the crossing predicted by Fig. [Hare exactly the same as those 
of actual levels of a maximally commuting operator with that ordering shown in Fig. [SJ The latter 
has been obtained by numerical diagonalization of a 5 x 5 maximally commuting operator (j2ip . 
H(u), with randomly chosen ^k, £fc and random dk obeying the above ordering. More examples 
of level diagrams are shown in Figs. El EJ and [9] and discussed in detail below. Next, we explore 
further consequences of Eq. (f46l) . 

Inevitability of level crossings for maximally commuting operators. In Sec. IIIII we have 
seen that the maximum allowed number of level crossings is {N — 1){N — 2)/2, see Eq. (j29|) . Now 
let us show that at least one crossing must be present in the spectrum of any maximally commuting 
H{u) = T + uV. Suppose the eigenvalues of V are ordered as di < dj < ■ ■ ■ < di < dm as shown 
in the diagram (j47p and assume there are no crossings. Then, the top level must go from i on the 
left to m on the right, i.e. i ^ m, the next level starting at j on the left must be connected to 
/ on the right, j — > I etc. Finally, we must have / -^ j and m ^ i. Consider in particular levels 
i ^ m and m ^ i. According to Eq. (|46p . this asymptotic behavior implies m = i — 1 (mod A^) 
and i = m — 1 (mod N). We obtain = 2 (mod N), which does not hold for any N > 3, i.e. the 
above assumption that levels do not cross cannot be true. Thus, we have demonstrated that at 
least one level crossing is always present. 

The maximum number of crossings M™^^ = {N — 1){N — 2)/2 is realized e.g. for the ordering 
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FIG. 5: Numerical energy levels of a 5 x 5 maximal Hamiltonian H{u) = T + uV with the same ordering of 



eigenvalues of V as that in Fig.H) Energies are scaled with a factor {t 



-I) 



-1/2 



to highlight their asymptotic 



approach to eigenvalues of V . Matrix elements of H{u) are generated using Eq. (PT|) with random 7^, Ek 
and random dk constrained to obey the ordering of Fig. ID Note that the number of crossings as well as the 
levels that cross are exactly the same as those predicted by Fig. H) 




FIG. 6: A schematic diagram corresponding to maximum level crossings with N — 5. As is evident, 
A^ -^ A^ — 1 line has N — 2 crossings, A^ — 1 — > A — 2 adds A — 3 new crossings and so on, till the line 3 — > 2 
adds only 1 new crossing. Thus, the maximum number of crossings, l + 2 + 3+... + (A^ — 2) = (A^— 1)(A — 2)/2, 
can be confirmed. 
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b) 






FIG. 7: All distinct level diagrams of iV = 4 maximal Hamiltonians, H{u) = T + uV, drawn according to 
(|T7)) . There are {N— 1)! = 6 distinct diagrams each corresponding to iV = 4 different orderings of eigenvalues 
of V (see the text). For example, diagram b) corresponds to di > da > ^4 > ^2 and three other orderings 
obtained with a shift of the indices by an integer mod N, e.g. d2 > d4 > di > d^ etc. The diagrams predict 
either one or three level crossings in agreement with Eq. (|45p and specify which levels cross. However, when 
the number of crossings is less than the maximum (three), additional multiple crossings of the same two 
levels can occur. This can increase the number of crossings by 2K, see Fig. [S] In the present case, the 
number of crossings for orderings b) and e) can increase from one to three. 



di > d2 > ■ ■ ■ > d^, see Fig. El In this case Eq. (j46|) implies N ^ N — 1 yielding m = N — 2 
crossings, A^ — 1 ^ A^ — 2 giving rise to another m = N — 3 crossings etc., so that altogether we 



have Y^ 



N-2 , 
m=l 



m 



level crossings. Note however that the schematic level diagrams, such as 



those shown in Figs. [H [6l [71 and[9l do not account for the possibility of two levels crossing more 
than once. For example, the level 2 ^ 1 in Fig. Uh can go below the level 1^4 and come back 
above it again generating two additional crossings, see Fig. [8l Therefore, multiple crossings of two 
levels can increase the total number of crossings Mc for a given ordering of dk by an even number 
except when Mc = M^^^. In the latter case, since Mc cannot exceed M™^^, multiple crossings 
of the same two levels are prohibited. We see that multiple crossings do not modify Eq. (|45p . 
Interestingly, numerically we have found that for as low as A^ = 8, multiple crossings of the same 
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0.5 rE 





b) 




FIG. 8: a) Numerical energy levels of a 4 x 4 maximal Hamiltonian H{u) = T + uV with the same ordering of 
dk as in Fig. [7)d. Energies are scaled with {u^ + 1)^^/^ to highlight their asymptotic approach to dk- Matrix 
elements of H{u) are generated using Eq. (|2ip with random 7^, £& and random dk constrained to obey the 
above ordering. The multiple (twofold) crossing of the top two levels increases the number of crossings from 
one as in Fig. [7)d to three. This illustrates the generic situation arising when multiple crossings increase 
the number of crossings by an even integer over and above the number enforced by the diagrams (|47p . 
Nevertheless, as discussed in the text, this does not affect Eq. (|i5)) . b) The schematic of a). 



a) 





FIG. 9: All distinct level diagrams for A^ = 3 maximal Hamiltonians, H{u) — T + uV, drawn according to 
(|47|) . There are ( A^ — 1 ) ! = 2 distinct diagrams each corresponding to A^ = 3 different orderings of eigenvalues 
of V (see the text). For example, diagram a) corresponds to di > d2 > rfs and two other orderings obtained 
with a shift of the indices by an integer mod A^, i.e. ^2 > da > di and d^ > di > d2. The diagrams 
predict a single level crossing in agreement with Eq. (|45p and specify which levels cross. Since this is also 
the maximum number of crossings for A^ = 3, no multiple crossings of the same two levels are allowed, cf. 
Fig. El and [i 
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two levels are very common. 

Thus far, we have established that the total number of crossings Mc in the spectrum of an 
arbitrary maximally commuting Hamiltonian is 1 < Mc < M™**^. By inspecting all level diagrams 
for 3 < A^ < 6, we have also found that for a given N the total number of level crossings changes 
in increments of 2 from M™^^ to 1 (2) for odd (even) M™^^, i.e. we verified Eq. (j45|) for these A''. 
Moreover, this equation is also supported by the parity considerations in the end of Sec. lIIII and is 
consistent with all preceding observations regarding the properties of level diagrams. As such, we 
adopt it without a formal proof. 

Let us also comment that cases when some of the eigenvalues of y, d^, are degenerate should 
be regarded as crossings at u ^ ±00. Equivalently, one can treat T and V on equal footing by 
defining H{u, v) = vT + uV . Then, degenerate dk correspond to crossings at v/u = 0, while the 
crossings considered above occur either at finite v/u or at v/u — > ±00, or equivalently at u/v = 0. 
For example, levels of the BCS Hamiltonian, which is a linear combination of Gaudin magnets[40t], 
Hbcs = ^ X]^«^*(-^) + const, cross at the value of the BCS coupling constant g = \/B = \/u = 00 



or, equivalently, at u = 0, see e.g. Refs 



.0 



and 



42. 



We conclude this section with a discussion of useful properties and examples of energy level 
diagrams. There are A^! diagrams for a given A^ corresponding to permutations of eigenvalues 
di, ^2) • • • ) dN- However, some of them are identical. Specifically, orderings di < dm < ■ ■ ■ < di and 
dj+a < dm+a < • • • < dij^a that differ by a shift of indices by an integer a yield identical diagrams, 
since Eq. (j46p is invariant with respect to the replacement k ^ k + a (mod A^). Because A^ different 
orderings can be generated using this shift, it leaves (A^ — 1)! distinct diagrams. For example, there 
are two distinct diagrams for A^ = 3, see Fig. [9l Each corresponds to three different orderings 
of dk- Both diagrams predict a single level crossing. Since this is also the maximum number of 
crossings for A^ = 3, repeated crossings of the same two levels are not allowed. Therefore, a single 
crossing of either two top or two bottom levels is the only option for A^ = 3 maximally commuting 
Hamiltonians. For A^ = 4 there are six distinct level diagrams shown in Fig. [71 Four of them - 
diagrams a), c), d) and f) in Fig. [7]" exhibit the maximum number, M™^^ = 3, of level crossings. 
In a manner similar to that of the A^ = 3 case this is the only option for the corresponding sixteen 
orderings of d^- In contrast, in diagrams b) and e) showing a single crossing, multiple crossings 
can occur. This will increase the total number of crossings from one to three, see Fig. [HI 
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VII. SUBMAXIMAL HAMILTONIANS 



The preceding sections have focused on maximally commuting Hamiltonians, where we have 
explicitly constructed these operators, solved them exactly, and used the solution to explain the 
level crossings in such systems. In this section, we explore Hamiltonians linear in a parameter u 
characterized by less than the maximum number of commuting partners. Most importantly, we 
demonstrate that some of these submaximal Hamiltonians have no energy level crossings, i.e. the 
inevitability of level crossings due to parameter-dependent commuting partners appears to be an 
exclusive property of maximal Hamiltonians. 

As discussed in the Introduction, a given family of maximal Hamiltonians contains A^ — 1 
nontrivial independent commuting operators (see the discussion in the paragraph following Eq. ([6])). 
It is reasonable to expect that there exist submaximal families with N — 2,N — 3 etc. Hamiltonians. 
Similar to Eq. ([5]), any common n- independent symmetry is assumed to be factored out by going to 
blocks of the same symmetry. We may adopt a convenient terminology, where families with N — L 
nontrivial Hamiltonians are identified as being Type L (cf. Type I and II of Ref. ly). Then, the 
maximally commuting Hamiltonians are Type 1, those with N — 2 commuting operators are Type 
2 and so on. Since a nontrivial family must contain at least two nontrivial commuting operators, 
the first nontrivial instance of Type 1 occurs for A = 3, Type 2 for A = 4 etc., where A is the 
dimensionality of the state space. 

First, let us construct 4x4 Type 2 Hamiltonians linear in a real parameter u. Our task 
is therefore to identify two 4x4 commuting real symmetric matrices that do not have the third 
independent commuting partner other than (a + ub)I. We will do so by employing the construction 
of maximal Hamiltonians detailed in Sec. |TI1 Consider 4x4 operators H {x,y,u) = xT + y K + uV 
and H {x,y,u) = xT + y K + uV, linear in parameters x, y and u, such that 



H{x,y,u),H{x,y,u) 



0. 



(48) 



Since this equation is to hold for all x, y, and u, the coefficients of the xy, xu, yu etc. terms must 
vanish individually. We obtain 



T,T 


= 


v,v 


= 


K,K 



0, 



(49) 



T,V 



T,V 



V.K 



T,K 



V.K 



T,K 



(50) 
(51) 
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FIG. 10: Numerical energy levels of a 4 x 4 submaximal Hamiltonian H{u) obtained from Eqs. ((52|) and 
([54| with xq — 1, random yo, dk, Efc, and Sk- Energies are scaled with {u^ + 1)^^^^ as in Fig. [H Unlike 
A^ = 4 maximal Hamiltonians, which always have two nontrivial commuting partners, this H{u) has only 
one such partner (see the text). Note that levels of H{u) do not cross at any u, i.e. the mere existence 
of a nontrivial commuting partner does not guarantee level crossings. This is to be contrasted with 4x4 
maximally commuting Hamiltonians which always exhibit either three or one crossings, see Figs. [7] and [H 



Let us choose these real symmetric matrices so that (T + uV, T + uV) and (K + uV, K + uV) 
belong to two different famihes of maximally commuting Hamiltonians parameterized by 7^ , £m 
and 7m = ^mi'^m-, respectively, see Eq. (J2T]) . Then, Eqs. (jl9]) and ([50]) are satisfied by construction. 
Consequently, it remains to solve Eq. (j5ip . 

Thus, from Eq. ()2ip we have in the common eigenbasis of V and V 
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where d\^ are the eigenvalues of V . Matrix elements of H {x,y,u) are obtained from Eq. (j52p by 
replacing d^ ^> d^. Using these expressions for the matrix elements, one can rewrite the remaining 
commutation relation (|5ip as follows 
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{Ei — Ej) {Ej — Ek) {Ei — Ek) {Ei — Ej) {Ej — Ek) {Ei — Ek) 

Therefore, choosing dr,dr,Er, and Er, we obtain 7^ from Eq. ([53|) . This yields two commuting 
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matrices H{x, y, u) and H{x, y, u). Fixing nonzero values of x = xq and y = yo, we obtain a Type 
2 family of Hamiltonians linear in u, 

H{u) = {xoT + yoK)+uV, H {u) = {x^f + y^K) + uV . (54) 

There are a number of equivalent ways to verify that these operators are indeed Type 2 rather 
than maximally commuting. For example, one can show that their matrix elements ()52p cannot 
be cast into the form ()2ip . Alternatively, it can be demonstrated that conditions (|A3p necessary 
for any maximal operator do not hold. However, a less formal, but more fruitful verification uses 
the following argument. We have seen in Sec. IVII that any N = 4 maximal Hamiltonian must have 
either one or three level crossings. Let us check if this holds for the Hamiltonians (|54|) . To this 
end, we set xq = 1, generate random yQ, dj., e^, and e^, and obtain 7^ from Eq. ([53]) and H{u) 
from Eq. (j52p . Doing so repeatedly and numerically diagonalizing the resulting Hamiltonians we 
observe that they always have either no or two level crossings. An example with no crossings is 
shown in Fig. [TOl Thus, operators ()54p are Type 2. 

We see that level crossings are not guaranteed when the number of commuting operators is 
less than the maximum - nontrivial solutions of Eq. Q do not necessarily imply crossings. The 
converse is also false, i.e. level crossings can occur in the absence of any nontrivial commuting 
partner linear in u and any u-independent symmetry. For example, one can show that the 4x4 
Hamiltonian given by Eq. ()A4p in the Appendix with a single level crossing at u = has no 
nontrivial commuting partners and no u-independent symmetry. Interestingly, A^ = 4 is the first 
dimensionality where this happens as for 3x3 real symmetric matrices linear in u a level crossing 
implies a nontrivial commuting partner linear in u and vice versa [18| . 

VIII. SUMMARY AND OPEN QUESTIONS 

In this paper, we addressed the problem of the violation of the Wigner-von Neumann non- 
crossing rule in quantum integrable systems. For this purpose, we introduced and studied a general 
class of maximal Hamiltonians - a vector space of N x N real symmetric Hamiltonians, H{u) = 
T + uV, characterized by the existence of the maximum possible number (N) of independent 
mutually commuting integrals similarly linear in the coupling u. We have resolved this commutation 
property and explicitly constructed general maximal Hamiltonians, see Eq. (j2ip . Interestingly, 
these operators are equivalent to the Gaudin magnets (j30p in the next to highest weight sector, 
J^ = J^i^^ — 1, where J^ is the z projection of the total spin. 
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The mapping to Gaudin magnets allowed us to obtain a complete exact solution for the eigen- 
states and eigenvalues of H{u), Eqs. (j40|) and (j41|) . Furthermore, we have demonstrated that energy 
level crossings are inevitable for maximal operators, i.e. there is always at least one crossing. The 
total number of crossings varies from 1 or 2 to (A^ — 1)(A^ — 2)/2, see Sec. IVII Thus, the mere 
existence of the maximum number of commuting partners guarantees a) an exact solution and b) 
level crossings. This relationship between the existence of conservation laws and exact solution 
is a quantum analog of the famous Liouville- Arnold theorem in classical mechanics. The latter 
states that if a classical model with n degrees of freedom has n Poisson-commuting integrals, its 



equations of motion are exactly solvable 2C ] . 



At the same time, by constructing an explicit example we have demonstrated that submaximal 
Hamiltonians - real symmetric operators of the form T + uV with less than the maximum number 
of linear in u commuting partners - often show no instances of level crossings at any u. Thus, 
the inevitability of crossings is an exclusive feature of maximally commuting operators. Similarly, 
we have also shown that the presence of level crossings does not necessarily imply the existence of 
a nontrivial commuting partner linear in n, i.e. crossings can occur even in the absence of such 
partners as well as n-independent symmetries. 

Our understanding of properties of parameter-dependent energy spectra in integrable models 
is far from complete. We conclude this section with a list of open questions stemming from the 
results of this work. 

1. We have shown that there are submaximal Hamiltonians with no level crossings. Never- 
theless, crossings often do occur in these systems in violation of the non-crossing rule, see 
Sec. lVIll This indicates that there is more to the link between crossings and the presence of 
commuting partners. It is interesting to understand this link for submaximal operators, what 
lifts level repulsion in this case, and why crossings happen only for a fraction of submaximal 
Hamiltonians. 

2. In Sec. IVIH we have also introduced a notion of Type L commuting family characterized by 
N — L nontrivial integrals. In this classification maximally commuting operators are Type 
1, while submaximal operators correspond to L > 2. A natural question is whether there 
is a general explicit parametrization for Type 2, 3 etc. similar to the one obtained in this 
paper for maximal systems. For instance, one can show that | J^| < <^max ~ ^ sectors provide 
examples of submaximal commuting families. 

3. What is the role of maximal Hamiltonians in the context of general quantum integrable 
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Hamiltonians? For the central spin Hamiltonians (Gaudin magnets) they represent the next 
to highest weight sector. Do other integrable models have maximally commuting sectors? 

4. In this paper, we focused on operators linear in the coupling u. An interesting question 
is how our results can be generalized to operators with a more general, e.g. polynomial, 
dependence on the coupling. 

5. We have established hard bounds for the number Mc of level crossings in A^ x A^ maximally 
commuting operators. Can one also determine the distribution of Mc for large A, i.e. the 
relative prevalence of maximal Hamiltonians with a particular number of level crossings? 
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APPENDIX A: SHASTRY'S CONSTRUCTION OF COMMUTING MATRICES 

In Ref . lis Shastry developed an approach to generate commuting real symmetric NxN matrices 
linear in a parameter u. Here we show that matrices obtained with this approach belong to 



maximally commuting set constructed in Sec. 
First, we briefly review the results of Ref. 



n 



see Eq. mh . 



19|. Consider Eq. ([3]). In the common eigenbasis of 



V and V the second relation in Eq. ^ becomes 

where Tij (Tij) are the matrix elements of T (T) and di and di are the eigenvalues of V and V, 
respectively. It remains to consider the [T,T] = commutation relation in Eq. ^. This can be 
cast into the following form: 

tJ-ijkdi + fJ-jkidj + fJ-kijdk + ^ J^iijkd-i = 0, distinct i,j, k, (A2) 

l^i,j,k 

where ^iijk and i^Hjk depend only on matrix elements of H{u) and not on those of H{u). Specifically, 
they involve only dr,Trr and Srm- 
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A set of particular solutions to Eq. ([3]) can be obtained by setting the coefficients in Eq. ()A2 
at each dr individually to zero, i.e. 



IJ-ijk = 0, uiijk = 0, distinct /, i,j, k, 



(A3) 



Now commuting H{u) and H{u) can be generated as follows. One first chooses 3A^ — 1 parameters, 
e.g. 2A^ — 3 variables Sir and S2r for {2,3} < r < N and A^ + 2 variables {d^}, Tn, and T22. Then, 
Eq. ()A3P reduce to linear equations and can be solved for the remaining variables. Once H{u) is 
determined in this way, H{u) can also be constructed, see Ref. |19i for details. This scheme is quite 
suitable for numerical implementation and, having examined several examples, Shastry observed 
crossings in all cases. Based on this and the results of Ref. [ISJ he conjectured that these matrices 
will always exhibit them. 

To show that this construction always yields maximal Hamiltonians, we note that Eq. ()A3|1 is a 
sufficient condition for Eq. ()A2p to have A^ linearly independent solutions for {di,d2, . . . disi). Since 
dr are the eigenvalues of V, the existence of A^ linearly independent solutions means that there are 
A^ Hamiltonians H{u) with linearly independent Vs. The absence of u- independent symmetries 
can also be demonstrated (it follows from Sij 7^ for all i 7^ j). Thus, we have a maximally 
commuting set, see the Introduction. The only difference is that by construction dr are not allowed 
to be degenerate, see Eq. (jAip . while the maximal set contains these matrices as well. 

Finally, we write down an example (see the discussion at the end of Sec. IVIB of a 4 x 4 
Hamiltonian H{u) with a level crossing but no u-independent symmetry and no commuting partners 
linear in u other than trivial ones - cH{u) + (o + ub)I, where a, 6, and c are real numbers and / 
is the identity matrix, 
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